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Abstract 

Within the one-dimensional tight-binding minibands and on-site Coloumbic 

interaction approximation, the absorption spectrum and coherent wavepacket 

time evolution in an optically excited semiconductor superlattice driven by 

dc-ac electric fields are investigated using the semiconductor Bloch equations. 

The dominating roles of the ratios of dc-Stark to external ac frequency, as well 

as ac-Stark to external ac frequency, is emphasized. If the former is an integer 

M, then also Af harmonics are present within one Stark frequency, while the 

fractional case leads to the formation of excitonic fractional ladders. The later 

ratio determines the size and profile of the wavepacket. In the absence of 

excitonic interaction it controls the maximum size wavepackets reach within 

one cycle, while the interaction produces a strong anisotropy and tends to 

palliate the dynamic wavepacket localization. 
PACS numbers: 78.47,+p, 73.50.Fq, 73.20.Dx 
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I. INTRODUCTION 



The effect of static and/or time-dependent electric fields on superlattices has been the 
subject of much recent work. Of particular interest is the possibility of terahertz emission 
caused by Bloch oscillations. HI] Focus in the past has been mainly on the presence of static 
electric fields while neglecting the THz emission signal propagating in the superlattices. 
Recently however, in order to interpret the experimental observation of peak shifts in the 
Wannier-Stark absorption lines of optically excited superlattices, a quasi-static model has 
been introduced where the generated THz electric field is treated adiabatically. [H Dignam 
has shown, however, that this adiabatic treatment of the THz electric field is inadequate. || 
The treatment can be regarded as roughly correct only with the proviso that the interband 
inverse dephasing time is greater than the THz field frequency. If this condition is not 
fulfilled it is then necessary to treat the static and time- dependent electric fields on an equal 
footing. |i 

It is also known that the combined dc and ac electric fields can induce many interesting 
and novel phenomena in superlattices. These include the appearance of inverse Bloch oscil- 
lators fHH, absolute negative conductivity , negative absorption (gain) 0, and fractional 
Wannier-Stark ladders [^]|| . In this paper we study the absorption and wavepacket dynamics 
in optically excited semiconductor superlattices under the influence of dc-ac electric fields. 
It is found that the parameters from both the dc and ac electric fields play important roles 
in determining their coherent dynamics behavior. The absorption spectrum is determined 
by the ratio of the Stark frequency ujb to the THz field frequency uo ac . When this ratio is 
an integer Af, M harmonics appear within the frequency interval ujb, while if it is a fraction 
of the form us/uJac = p/q the number of harmonics is q, each being identified as part of 
the so called fractional ladder. It is the first time that this fractional structure is found in 
the frequency domain by taking into account the excitonic interaction. Also, the excitonic 
together with Fano interference on the asymmetric lineshape of the absorption is pointed 
out. On the other hand, the extent and profile of an electron- hole wavepacket is mainly 
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determined by the ratio of the Stark frequency associated with the strength of the THz 
field, and its frequency. The excitonic interaction produces anisotropy in the wavepacket 
profile and reduces the dynamic localization effect. 



II. MODEL AND NUMERICAL RESULTS 

The coherent dynamics of an optically excited semiconductor superlattice (SL) driven 
by static and time-dependent electric fields can be described by the following semiconductor 
Bloch equations 



^ + ~F(t) • Vk) Pu(t) = -^Kk + e^k - iV L ]P k (t) - ~[n e>k + n h>k - l}u R , k (1) 



^ ± |F(t) ■ Vk^) n e(h)k (t) = -2Im[oj R)k P*] - T T n e , (h)k (t) , 
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where P k (t) is the interband polarization and n e ^ h ) k (t) the electron (hole) population density 
in the conduction (valence) band. The quantities e i)k = e i k — S q ^|k-qpj,q, (i = e, h) are 
the renormalized electron and hole energies due to the Coulomb interaction, with e^k being 
the single-particle band energy. Also, ujR tk = d cv £ + J2k M k -q|-^q are ^ ne ^normalized Rabi 
frequencies, with d cv the dipole moment and £{t) = £ exp(—t 2 /r 2 ) exp(-iuLt) (r = 100 fs), 
the Gaussian laser pulse profile. The relaxation time approximation has been assumed, with 

and IV being the longitudinal and transverse relaxation rates, respectively. 

We want to obtain the spectrally-resolved absorption coefficient and coherent wavepacket 
dynamics for an optically excited SL driven by the combined dc-ac fields F(t) = Fq + 
Fi cos(cu ac t). The both components of the electric field are assumed to point along the same 
direction. 

For simplicity, a one- dimensional tight-binding (TB) approximation has been adopted 
with the two miniband widths and gap energy being A c (conduction), (valence) and E g 
respectively. [f^ |TT| , |i"3"| This approximation neglects the in-plane continuum contribution to 



the absorption, whose consequence is discussed in the following. Firstly, in the absence of 



excitonic interaction, the absorption a{uj) is proportional to the 3-D joint density p( 3 )(co>), 
which is the convolution product of the one dimensional joint density of states (JDOS) p^(e) 
and 2-D continuum JDOS p^(e). The latter is a step function changing abruptly at e = 0. 

a(u) oc / p (2) (e)p (1) (^-e) de. 



The one dimensional JDOS p^(e) looks like a trapezoid with the concave part in the middle, 
and it approximately extends between [—A, A], where A here is the combined miniband 
widths defined as A = Ac + A^. Due to delayed kernel of p^'{fuo — e) in the above expression, 
the resulting p^(cu) begins to attain the finite value at u — —A/h. and becomes an 
approximate constant at uj = A/h and beyond. That is to say, that the behavior of p^(e) at 
e = changes completely due to one-dimensional JDOS. Thus, the inclusion of the in-plane 
continuum contribution can influence the height of characteristic peaks. Because JDOS is 
always positive, it can not turn the gain-type lineshape to the resonance-type shape or vice 
versa. Secondly, the contribution to the oscillator strength from the in-plane continuum 
is negligible in the range of terahertz radiation of experimental interest jTjynj. Thirdly, 
one also can regard the in-plane continuum as the source of the inhomogeneous broadening 
which was proposed by Thomas and von Plessen [Hjj . This inhomegeneous broadening may 



smear some characteristics peaks appearing in our present numerical simulations. 

This TB approximation captures the qualitative features of the realistic superlattices, 
and can give a fairly good description of their response within the strong coupling regime. 
|T3| It has been used with success by many groups in the past. For the Coulomb 



interaction we adopt the on-site model. In the absence of the external electric fields, the 



contact Coloumb potential produces single bound states [15,18]. In the dc -£LC CclSC, however, 



the situation changes significantly. For example, in the realm of the near resonance case 



of stark and ac field frequencies, unequal spacing |19| between the Wannier-Stark ladders 
due to excitonic interaction can make the quasienergy of the system to present more than 
one bound Floquet state. The number of bound Floquet states is determined by how many 
distinct deviations between excitonic Wannier-Stark ladders and multiple ac field quantum. 



This kind of Floquet bound states usually appears in the low Wannier-Stark ladder indices. 
This case of more than one bound Floquet states has been discussed in great detail by Liu 
and Zhu clearly within the picture of dressed Wannier-Stark ladders [p|,|T5|j. 



Because we assume a weak ultrashort laser pulse to excite the SL, instead of the pertur- 
bative expansion of semiconductor Bloch equations in terms of the order of So, we take into 
account the full contributions from interband polarization. The nonlinear effect facilitated 
by the coupling of n(k,t) with p(k,t) on the lineshape of the fractional ladders will be dis- 
cussed in the following text. It is well known that the absorption coefficient can be expressed 
as a{uj) oc lmx{uj) = lm[P(u)/ '£(ui)], with P(u>) and S{oS) being the Fourier-transformed 
quantities of P(t) = J2kPk(t) and S(t), respectively. We discuss separately two cases of 
interest, in which the ratio of the Stark frequency ujb = eF^d/h to the THz field frequency 
ujac is (a) an integer TV, and (b), an irreducible fraction p/q. 

Equations (1) and (2) were numerically integrated in the accelerated k basis, where the 
drift term can be eliminated. Our results for the integer case ujb = 2uj ac = 7ir THz are 
shown in Fig.l. Figure 1(a) is for the static case, while Fig.l(b) is for eFid/huo ac = 0.1 . 
The values of parameters used are: combined miniband width A = A c + A„ = 40 meV; 
Coulomb interaction strength V =10 meV, comparable to miniband width; dephasing rates 
r-r = = 0.5 TH; In the simulation of the coherent time evolution of the wavepackets, 
the central frequency of the ultrashort laser pulse ujl is tuned to be in resonance with the 
energy-gap frequency u g . 

The Wannier-Stark resonances are clearly identified in Fig. 1(a). The frequency spac- 
ings are not equal, however, due to the excitonic interaction [|19j. In Fig. 1(b) the spec- 



trum appears more complex, showing a weak additional structure between any two main 
peaks. These weak structures exhibit either gain or absorption lineshapes, resulting from 
the nonlinear sum and difference frequency mixing of the dc and THz electric fields. It 
should be emphasized here that the ratio edF\j7iio ac is a dominating factor in determin- 
ing whether gain or absorption occurs. That is to say, changing the ratio edFi/%uo ac may 
turn the resonance-type peaks into gain-type dips, or vice versa. The linewidth broaden- 



ing of these structures also depends on the ratio edFi/Twj ac . This can be understood by 
noting that the quasi-energy band [^,^1|] of a dc-ac driven single-band TB system [F22| is 
e(k,u ac ) = + ^J n (edFi/huj ac )cos(kd), with ujB/oj ac = n. According to the general- 
ized golden rule in the single-particle picture the transition rate P c ^ v (k) from the valence 
dressed band to the conduction dressed band for every component k, is proportional to 
J2 n 3[ e c(k, oj ac ) — e v {k,uj ac ) + nhu) ac — Tuol)] • HU For every k, the transition can be viewed 



as a two-level system (TLS) with different energy spacing. These different energy spacings 
together with the k-mixing effect of the excitonic interaction lead to the broadening of the 
resonance peaks and gain dips. Of course, one can select the ratio edFi/Tiuj ac to control the 
quasi-energy band widths, hence the broadening. Another feature, apparent when compar- 
ing Figs. 1(a) and 1(b), is the relative blueshift of the main resonance peaks in the latter. 
These shifts are not equal, as emphasized by the horizontal arrows. Other integer cases such 
as Af us/uJac = 3, 4 were also calculated, yielding Af resonance-type peaks or gain-type 
splittings in every Stark frequency interval wg. In order to save space, they will not be 
discussed here. 

Distinct asymmetric lineshapes of the peaks appear both in dc/dc-ac field case in Fig.l. 



This asymmetry can be attributed to the Fano resonance |24j|. The theoretical prediction of 



the Fano resonance in dc-biased semiconductor superlattices was made by Whittaker ||25|| . 



Linder [IS], Glutsch, and Bechstedt [26], and was experimentally observed by Holfeld et al 



1 27J • The mechanism of the Fano resonance is due to the coupling between the discrete and 
continuous excitonic states of the relative motion, mediated by the Coulomb interactions 
25| , |26|jl9| , |27[| . The Fano resonance in dc-ac (ujb/uj an integer) driving superlattices were 



studied by Liu and Zhu |15|]. The resonance takes effect through the coupling between the 
discrete quasienergy excitons and sidebands of their continua. 

The other interesting case is when the ratio ujB/u ac equals an irreducible fraction p/q. 
Figure 2 shows the spectrally- resolved absorption coefficient for the cases UB/uj ac = 1/3, 2/3 
and 4/3. Here the same excitonic interaction strength and combined miniband width as in 
Fig.l are used. In all cases edF\/Tiuj ac = 1.0, and uo ac is fixed to be 7n THz. 



Surprisingly, the resonances appear at approximately the same frequencies in all three 
panels, although the Stark frequency in the top and bottom panels differ by a factor of four. 
In order to guide the eyes, dash-dot-dot vertical lines were drawn around these resonances. 
It can also be seen from the figure that some resonances split into two peaks, which can be 
understood in the following context. The beating frequency between the peaks in all panels 
is almost the same: |o; oc . This frequency equals neither the Stark frequency ujb, which 
changes from panel to panel, nor the angular frequency of the ac field u ac . As a matter of 
fact, the quasi-energy band of the Hamiltonian driven by dc-ac electric fields with the ratio 
^s/^ac = p/q, is the underlying mechanism for the appearance of these peaks. Within the 
one-electron and single-band approximations, these quasi-energy bands are given by 

Tl 

e n (k) = A q (k,u B ) + -haj B , n = 0, 1, • • -,p - 1 , (3) 

where A q (k, Ub) gives the band dispersion. || From the last term in the right hand side 
of Eq. (3) it is evident that the quasi-energy band consists of p subbands with equal energy 
spacing hu B /p(= hu ac /q). The absorption spectrum reflects this energy gap through the 
generalized Fermi golden rule described above, while the term A q (k,u B ) is responsible for 
the broadening of the spectrum. The different values of u B taken in the three panels result 
in different values of A q (k,UB), hence, different broadening. Also, the inclusion of the 
Coulomb interaction makes the positions of the peaks to not exactly match the w ac /3 rule. 
Yet, since on the average there is one peak per u> ac /3 frequency interval, the structures are 
termed excitonic fractional ladders. Their temporal counterpart has been discussed recently 



in the time-resolved four- wave mixing spectrum ||17|| . Results in the absence of excitonic 
interaction are presented in Fig. 3, where we used the same parameters as in Fig. 2 except 
that now V = 0. The fractional Stark ladder of spacing u) ac /?> can be clearly identified, 
with peaks distributed symmetrically about uj = u g , matching the co> ac /3 rule better than 
the excitonic case. These have been indicated by the vertical arrows in the upper panel 
of Fig. 2. Comparison of the fractional ladder lineshape in Fig. 2 and Fig.3 reveals that the 
purely Lorentzian lineshape of some peaks in Fig.3 is absent in the corresponding peaks 



in Fig.2, where the excitonic interaction is taken into account. This demonstrates that 
the excitonic interaction together with the Fano resonances destroys the symmetry and 
smears some resonance peaks. Here, the Fano interference with small parameter qj takes 
effect through the coupling of the discrete excitonic fractional ladders with the continuous 
sideband. It should be emphasized that the absorption signal depends also on the phase 
difference between the ultrashort laser pulse and the ac field. This important effect has been 
found recently by Meier [28| et al in an anisotropic 3-D superlattice. 



When the exciting laser field becomes strong, the nonlinear effect due to the high order 
exciting laser field magnitude 8 begins to take effect. The typical simulation of the effect is 
through the pump-probe configurations, where the weak probe laser field is delayed to time 
t p with respect to the strong pump field. Then the different linear and nonlinear signals can 
be extracted [[29|,[3(J]. This nonlinear effect makes the clear-cut fractional ladder harmonics 
broadened, tending to merge with each other. This is ahown in Fig. 4, where the ratio of the 
exciting laser amplitude has been indicated in the upper-right corner of the respective panels. 
Here we only show the case of uos/^ac = 1/3, and we use otherwise the same parameters as 
those in Fig. 3, setting the time delay r p to 50 femtoseconds. 

There is another interesting but subtle problem when the ratio UB/uj ac is an irrational. 
Numerical calculation can not identify the essential difference between the rational and 
irrational cases. Physically, this difference will also be blurred by the scattering processes. 

We now turn to the coherent dynamic evolution of electron- hole wavepackets (WP). 
Results for a pure dc field have been previously reported by Dignam and coworkers, where 
Bloch oscillations are clearly demonstrated. |14[ Meier et al studied relative motion in the 



exciton wavefunction in momentum space and found a dimensional crossover due to dynamic 



localization induced by a strong THz field. [pl],[33] More recently Hughes and Citrin gave 



vivid pictures of the strongly anisotropic excitonic WP in optically excited semiconductor 
quantum wells driven by an in-plane THz field, through a time sequence of the WP snapshots. 



The WP coherent evolution may be described through the squared modulus of the real 



space Fourier transform of the interband polarization P(z,t) = | J2k ex P( — ' z )P(k, t)\ 2 , 
giving the real space description of the probability density of finding an electron and hole 
separated by z at time t. [|33| It should be noted that wavepackets can not directly be probed 
by experiments, although it is theoretically interesting. In optically excited semiconductors, 
the center of mass quasi- momentum of the electron- hole is always zero || , which makes the 
description of the behavior of the wavepackets in the center of mass difficult. 

In order to elucidate the role the excitonic interaction plays on the profile of wavepackets, 
we first present results in the absence of this interaction. Snapshots at various times are 
shown in Fig. 5, where the Stark and ac frequencies are chosen to have the same value 4n 
THz (T ac = 500 fs), and eFid/?UJ ac is set close to the first root of the first order ordinary 
Bessel function (3.832). 

Note that at all times the WP is symmetric about the position z = 0. With the elapse of 
time, it first splits into two parts and subsequently into three parts at half T ac , when the size 
of the WP acquires its largest value within a cycle. After T ac /2, the size begins to reduce. 
Due to the synchronous dynamic localization of electron and hole, at time t = T ac , the size 
of the WP experiences the severe shrinkage shown in the last panel of Fig. 5. At later times, 
the WP will undergo the same expanding and shrinking processes in every time cycle (T ac 
). This phenomenon demonstrates that dynamic localization has the ability to control the 
size of the WP. That is to say, that dynamic localization will make it impossible to find the 
electrons and holes simultaneously beyond a certain value of the relative coordinate. Also, 
comparison Fig. 5 and Fig. 6 shows that dynamic localization makes the value of P(z,t) 
larger at z = 0, demonstrating that it will enhance absorption. 

Figure 6 illustrates the same case as that of Fig. 5, except that we choose eFid/hu ac to 
be 4.5. Now the WP also splits into three parts at time t = T ac /2. The most important 
difference, however, lies in that now the size of the WP does not experience shrinking 
towards the end of the cycle. In fact, the size of the WP grows steadily with time, until 
it disappears due to dissipation. We next turn the excitonic interaction on, with results 
showing a substantial change in the profile of the WP, as shown in Fig. 7. Here we use 
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the same parameters as in Fig. 5. The WP becomes strongly anisotropic. Its height is 
reduced and although its extent still undergoes a shrinkage towards the end of the cycle, this 
shrinkage is less severe as compared with the non interacting case. The splitting structure 
can be viewed as the spatial interference |33[] induced by the coulombic re-scattering due to 
the oscillating motion of the photo-generated electron-hole driven by dc-ac fields. 



III. CONCLUSIONS 

In summary, the spectrally-resolved absorption and real-time wavepackets of an optically 
excited semiconductor superlattice driven by combined parallel dc-ac electric fields are ob- 
tained by numerical integration of the semiconductor Bloch equations. The ratio of dc field 
Stark frequency ujb to ac field frequency u ac determines the characteristic peaks of the ab- 
sorption spectrum. When the ratio is a fraction, we identify for the first time the fractional 
ladders in the energy domain by taking into account the excitonic interaction, which we take 
as comparable to the miniband width. For the integer M case we find M features within 
the frequency interval ub, corresponding to gain and resonant absorption. The effect of the 
dc-ac field on the size and profile of the WP shows that, as expected, dynamic localization 
vigorously curtails the elongation of the WP within one cycle makes its evolution behavior 
to be of breathing mode. The excitonic interaction destroys the WP symmetry about its 
center and softens the dynamic localization of the WP within one time cycle. 
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FIGURES 

FIG. 1. The spectra-resolved absorption coefficient a(uo) for uos/^ac = 2 is plotted as a 
function of to — ui g for two cases (a) edFi/fioj ac = 0.0 (pure dc); (b) edFi/fiuj ac = 0.1. 
the additional harmonics appears between two neighboring Stark resonances. 

FIG. 2. The spectra-resolved a(oj) for the fractional cases uos/^ac = 1/3,2/3,4/3 is plotted as 
a function of u> — oo g , showing the excitonic fractional ladders. The ratio edFi/Tiu) ac is chosen to be 
1.0, while the other parameters are the same as in Fig.l 

FIG. 3. The same as Fig. 2, except that the excitonic interaction is absent. The fractional 
ladders are now symmetric about u g 

FIG. 4. The pump-probe spectrally-resolved four-wave mixing signal squared |P(w)| 2 has been 
shown in the four panels with the pump laser field magnitude ratio shown in the right-upper corner. 
Only the case of uj ac = uj 9 =1/3 is shown, the other parameters used here is the same as those in 
the Fig.2. 

FIG. 5. Coherent wave packet evolution in time in the absence of excitonic interaction. The 
ratio edFi/htOac is chosen to be 3.832, leading to dynamic localization within one cycle. 

FIG. 6. The same as Fig.4, except that the ratio edFi/Uu) ac is chosen to be 4.50, leading to 
absence of dynamic localization 

FIG. 7. The same as Fig.4, except that excitonic interaction is taken into account 
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